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Abstract 

The Bethe-Salpeter equation restores exact elastic unitarity in the s— channel 
by summing up an infinite set of chiral loops. We use this equation to show how a 
chiral expansion can be undertaken by successive approximations to the potential 
which should be iterated. Renormalizability of the amplitudes in a broad sense 
can be achieved by allowing for an infinite set of counter-terms as it is the case in 
ordinary Chiral Perturbation Theory. Within this framework we calculate the vrvr 
scattering amplitudes both for s— and waves at lowest order in the proposed 
expansion where a successful description of the low-lying resonances {a and p) and 
threshold parameters is obtained. We also extract the SU (2) low energy parameters 
^1,2,3,4 from our amplitudes. 
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Figure 1: Diagrammatic representation of the BSE equation. It is also sketched the used kinematics. 

1 Introduction 



Chiral Perturbation Theory (ChPT) to finite order is unable to describe resonances. Actu- 
ally, it is rather the opposite, resonances determine the bulk of the 0{p'^) parameters 0]. 
Such a description requires the use of a non-perturbative scheme. Several approaches have 
been suggested. Fade Re-summation (FR) 0, Large iV/— Expansion (LNE) |Q, Inverse 
Amplitude Method (lAM) P], Current Algebra Unitarization (CAU) P], Dispersion Re- 
lations (DF) 1^, Roy Equations Coupled Channel Lippmann-Schwinger Approach 
(CCLS) and hybrid approaches |10[- Besides their advantages and success to describe 
the data in the low-lying resonance region, any of them has specific drawbacks. In all 
above approaches except by LNE and CCLS it is not clear which is the ChFT series of 
diagrams which has been summed up. This is not the case for the CCLS approach, but 
there a three momentum cut-off is introduced, hence breaking translational Lorentz in- 
variance and therefore the scattering amplitude can be only evaluated in the the Center 
of Mass (CM) frame. On the other hand, though the LNE and CAU approaches preserve 
crossing symmetry, both of them violate unitarity. Likewise, those approaches which 
preserve exact unitarity violate crossing symmetry. 

In this paper we propose the use of the Bethe-Salpeter equation (BSE) to sum up 
an infinite set of diagrams without use of cut-offs in the physical amplitudes. We will 
show that this re-summation restores exact elastic unitarity in the s— channel and it 
naturally leads to the appearance of the experimentally observed resonances. Besides, 
crossing symmetry can be restored perturbatively. Our approach, at lowest order and in 
the chiral limit reproduces the bubble re-summation undertaken in Ref. |]1T 



2 The Bethe-Salpeter Equation. 

For the sake of simplicity, we neglect coupled channels contributions and thus the BSE for 
the scattering of two identical pseudo-scalar mesons of mass m and kinematics described 
in Fig. |1|, reads 

Tp{p,k) = Vp{p,k) + i [ -^Jp{q,k)A{q + P/2,m)A{-q + P/2,m)Vp{p,q) (1) 
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where Tp{p,k) and Vp(p, fc) are the total scattering amphtudeQ and the two particle 
irreducible amplitude respectively. Besides, A is the exact pseudoscalar meson propagator. 
Note that, to solve the above equation both the off-shell potential and amplitude are 
required. Clearly, for the exact potential V the BSE provides an exact solution of the 
scattering amplitude T |]T^. Obviously an exact solution for T is not accessible, since 
V and A are not exactly known. We propose an expansion along the lines of ChPT 
both for the exact potential {V) and the exact propagator (A). Thus at lowest order in 
this expansion, V should be replaced by the 0{p'^) chiral amplitude ('^^-'T) and A by the 
free meson propagator, A°(r, m) = (r^ — + ie)~^. Even at lowest order, by solving 
Eq. (|l]) we sum up an infinite set of diagrams^]. Our expansion is related to the approach 
recently pursued for low energy A^A^— scattering where higher order t— and u— channel 



contributions to the potential are suppressed in the heavy nucleon mass limit [O]. 



To illustrate the procedure, let us consider elastic tttt scattering in the s— and p— waves. 
There, for comparison with the experimental phase shifts, 5/j(s), we define the projection 
over each partial wave J for each isospin channel I 

Tu{s) = l rPj{cose)TUp,k)d{cose) = '^""l ,, [e'''"^'^-l] (2) 

— * 

where 9 is the angle between p and k, Pj are Legendre polynomials and A(a;, z) = 
+ + — 2xy — 2xz — 2yz. 



2.1 Isoscalar s— wave tttt Scattering. 

At lowest order, the off-shell potential V in this channel is approximated by 

V^{p,k)^(^^n{P,k) = + (3) 

where m is the pion mass, for which we take 134.98 MeV, and / the pion decay constant, 
for which we take 92.4 MeV. To solve Eq. (0) with the above potential we propose a 
solution of the form 

T°(p,A;) = A{s) + B{s){p' + e) + C{s)p^e (4) 

where A, B and C are functions to be determined. Note that, as a simple one loop 
calculation shows, there appears a new off-shell dependence {p^k"^) not present in the 
0{p^) potential ^'^^T. That is similar to what happens in standard ChPT 

-'^The normalization of the amplitude T is determined by its relation with the differential cross section 
in the CM system of the two identical mesons and it is given by dcr/dil = \Tp{p, fc)p/647r^s, where s = P^. 
The phase of the amplitude T is such that the optical theorem reads lmTp{p,p) — — (itotls^ — 4sto^)^/^, 
with (Jtot the total cross section. 

^Note that our approach is different to the LNE one There the bubbles are summed up by means 
of a collective field with no well defined isospin. Hence, crossing symmetry is maintained but elastic 
unitarity in the s— channel is violated. 



2 



The above ansatz reduces the BSE to a hnear algebraic system of equations which 
provides the full off-shell scattering amplitude. The resulting inverse amplitude on the 
mass shell and in the CM frame (P = 0,p^ = = 0, P° = ^/s) reads 

^-1.^ ^ , 2(/^ + /,(4m^)f 

°° ^ ' °^ ^ 2/4(4m2) + (m2 - 2s) /2 + (s - 4m2)/2(4m2) ^ ^ 

where 

J(.-[J^ (Q. 

'"^'^ ' J (27r)4 [(^ _ £)2 _ ^2 + i,] [(g + £)2 _ ^2 + i,] ^ ^ 

/o,/2 and J4 are ultraviolet divergent integrals and thus Eq. requires renormalization. 
We have made use of translational and Lorentz invariance which relate the integrals his) 
and liis) with /o(s) and the divergent constants J2(4m2) and I^iAm'^). Note also that 
Jo(s) is only logarithmically divergent and it only requires one subtraction, ie., /o(s) = 
/o(s) — /o(4m2) is finite and it is given by 



1 / 4m2, + 1 



where the complex phase of the argument of the log is taken in the interval [— vr, 7t[ 

To renormalize the amplitude given in Eq. we note that in the spirit of an Ef- 
fective Field Theory (EFT) all possible counter-terms should be considered. This can be 
achieved in our case in a perturbative manner, making use of the formal expansion of the 
bare amplitude T = V + VGqV + VGqVGqV + ■ ■ ■, where Gq is the two particle propaga- 
tor. Thus, a counter-term series should be added to the bare amplitude such that the sum 
of both becomes finite. At each order in the perturbative expansion, the divergent part 
of the counter-term series is completely determined. However, the finite piece remains ar- 
bitrary as long as the used potential V and the pion propagator are approximated rather 
than being the exact ones. Our renormalization scheme is such that the renormalized 
amplitude can be cast, again, as in Eq. (|]). This amounts in practice, to interpret the 
previously divergent quantities /o,2,4(4m2) renormalized free parameters. After having 
renormalized, we add a superscript R to differentiate between the previously divergent, 
-^0,2,4(4^^), and now finite quantities, J(f2 4(4m2). These parameters and therefore the 
renormalized amplitude can be expressed in terms of physical (measurable) magnitudes. 
In principle, these quantities should be understood in terms of the underlying QCD dy- 
namics, but in practice it seems more convenient so far to fit /(f2 4(4m^) to the available 
data. The threshold properties of the amplitude (scattering length, effective range, etc..) 
can then be determined from them. Besides the pion properties m and /, at this order 
in the expansion we have three parameters. The appearance of three new parameters is 
not surprising because the highest divergence we find is quartic (his)) and therefore to 
make it convergent we need to perform three subtractions. This situation is similar to 
what happens in standard ChPT where one needs to include at the one loop level some 
new low-energy parameters (k) |Jl|]. In fact, if t— and m— channel unitarity corrections are 
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Taylor expanded around threshold, a comparison of our (now) finite amplitude, Eq. (^), 
to the O^p"^) mr amplitude in terms of these parameters, ^i, 2,3,4, becomes possible. We do 
this explicitly in the Appendix. 

At the lowest order in the expansion proposed in this work, we approximate, in the 
scattering region s > 4m^, the (9(1//^) t— and u— channel unitarity corrections (function 
hij in Eq. (|15|)) by a Taylor expansion around threshold to order (s — 4m^)^. At next 
order in our expansion (when the full (9(p^)-corrections are included both in the potential 
and in the pion propagator) we will recover the full t— and u— channel unitarity logs at 
0(1//^), and at the next order {0{l/f^)) we will be approximating these logs by a Taylor 
expansion to order (s— 4m^)^. Thus the analytical structure of the amplitude derived from 
the left hand cut is only recovered perturbatively. This is in common to other approaches 
(PR, DP, lAM ■ ■ ■) fullfiling exact unitarity in the s-channel, as discussed in [B, |0],[14 . 



2.2 Isovector j9— wave tttt Scattering 

At lowest order, the off-shell potential V in this channel is approximated by 

vxp,k)^^'^n{p,k) = ^ (8) 

As before, to solve Eq. (|l]) with the above potential we propose a solution of the form 

T^{p,k) = M{s)p - k + N{s){p- P){k- P) (9) 

where M and are functions to be determined. Note that, as expected from our previous 
discussion for the s— wave, there appears a new off-shell dependence {{p ■ P){k ■ P)) not 
present in the O(p^) potential. Again, this ansatz reduces the BSE to a linear algebraic 
system of equations which provides the full off-shell scattering amplitude. The resulting 
inverse amplitude on the mass shell, after angular momentum projection, and in CM 
frame reads 



T^^ (s) = = -Io(s) H 

s — 4m^ 



11 = = -^o{^) ^ — : — 7—^ — (10) 



Similarly to the s— wave case, the above equation presents divergences which need to be 
consistently removed in terms, for instance, of the scattering volume and effective range 
in the p-wave. Because the highest divergence present is quadratic only two subtractions 
are needed and hence only two undetermined parameters appear. As we will discuss 
below, crossing symmetry at (9(1//^) provides a relationship between the parameters 
entering in different isospin channels. We emphasize once more, that the proliferation of 
undetermined parameters is not a drawback as compared to standard ChPT. 
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Figure 2: s— (left) and p— wave (right) tttt phase shifts as a function of the total CM energy y/s. Left 



panel: Triangles, squares and crosses stand for the experimental analysis of Refs. |16{ and |17 



respectively. Right panel: Triangles and crosses stand for the experimental analysis of Refs. ]18| and 1 19 
respectively. Solid lines are the fit (see Table |l|) of our s-channel unitarized ©(p^)— model to the data of 
Refs. and lil. 



3 Results 



The solutions of Eqs. (^ and (|T0|) satisfy elastic unitarity. Indeed, the imaginary part of 
Tfj{s) is determined by the imaginary part of Io{s) and in the physical region is given by 

ImT7/(s + ie) = -ImJo(s + ie) = \' (11) 

in agreement with the unitarity requirement provided m is the physical pion mass. That 
makes possible to extract the phase shifts unambiguously within our framework. 

As we have already mentioned, our approach violates crossing symmetry. At order 
0(1//^) our isoscalar s— and isovector p— wave amplitudes are polynomials of degree two 
in the variable (s — 4m^), with a total of five (3+2) arbitrary coefficients (see Eq. (|T7p), 
and there are no logarithmical corrections to account for t— and m— channel unitarity 
corrections [hu term in Eq. ([T3|)). Far from the left hand cut, these latter corrections can 
be expanded in a Taylor series to order (s — 4m^)^, but in that case the one loop SU{2) 
ChPT s— and isovector p— wave amplitudes can be cast as second order polynomials in 
the variable (s — 4m^), with a total of four (/i,2,3,4) arbitrary coefficients |1[]. To restore, 
in this approximation, crossing symmetry in our amplitudes requires the existence of 
a relationship between our five undetermined parameters. This relation reads (see the 
Appendix) 

75/2^'72m2 + 8/o^'^ + 33/o'''^ + 5/f Vm^ + = (12) 

Once this constraint is implemented in our model, there exists a linear relation between our 
remaining four undetermined parameters and the most commonly used /i 2 3 4 parameters 
(Eq. (pp). 
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In Fig. H we show the agreement of our model with the experimental phase shifts, both 
in the s— (left) and (right) waves. We fit the four undetermined parameter^ I^'l^Am?) 
and /((^2'(4m^) to the scalar and vector data. Results of the combined fit can be found 
in Table |l]. Values of x^/dof and the threshold parameters deduced from our formulae 
are also shown in Table |l]. As we see, the vector channel is satisfactorily described up to 
1 GeV, whereas the scalar channel is well reproduced up to 0.8 GeV. In the latter case, 
and for these high energies, one should also include the mixing with the KK channel as 
pointed out recently in Refs. l^-ll^. Regarding the deduced threshold parameters we 
find agreement within experimental uncertainties. For completeness we quote the deduced 
values for the low energy coefficients (see Eq. ( [T^ ) in the Appendix) 

= -0.90 ± 0.09, k = 6.03 ± 0.08, = 2 ± 7, k = 3.9 ± 0.4 (13) 



if a fit to the data of Refs. |T^ , [|I^ is performed^. As we see they are in reasonable good 



agreement with those obtained in standard ChPT which also use the same set of data 
to analyze vrvr scattering. When the data of Ref.||l9l for the p— wave are fitted we find 



compatible values for /1.2 and incompatible values for ^3^4. We prefer to use the data of 
Ref. ||18| to estimate the low energy parameters ^i, 2,3,4? since they come much closer to 
threshold than those of Ref. ||T^ {^/s = 350 MeV versus 450 MeV). Furthermore both 



sets of data are incompatible up to about = 650 MeV. 

Once we have determined the parameters /i, 2,3,4, the / = 2, J = lowest order BSE 
amplitude is completely fixed, as demanded from crossing symmetry. It is to say, the 
corresponding integrals which appear in this channel are completely determined by the / 's 
or equivalently by the a— and p— channel integrals. We have evaluated the corresponding 
phase-shift in the I = 2, J = channel, using the values of Eq. (0), and compared to the 
data of Ref. |20| finding a remarkable agreement up to about 1450 MeV (x^/ (num. data) = 
0.9) and threshold parameters (maao = -0.0390±0.0018, and mHao = -0.0701 ±0.0010) 
in agreement within uncertainties with the experimental values. 



Comparison to recent two-loop calculations ||21| would also be possible. To make such a 
comparison really meaningful, we should include in our analysis higher waves than the s— 
and p— waves considered here, which would necessarily require to include 0{p^) corrections 
to both the potential and the pion propagator. This point is out of the scope of this work 



and will be subject of future research p2[. Nevertheless, we would like to point out that 
the lowest order explored in this letter, though very simple, it already leads to reasonable 
sizes for the two loop contributions to the amplitudes. Thus for instance, for the s— wave 
scattering length we obtain aoo/oocf*' = 1 + 0.27 (one loop) + 0.07 (two loops) + ■ ■ ■, in a 
reasonable agreement with the findings of Bijnens et al, in Ref. [pT[|. 



Recent work, Refs. [§]-[|T0]5 suggests that off-shellness can be ignored when solving 
the BSE, since it only amounts to a pion mass and pion decay constant renormalization. 
This on-shell procedure corresponds to setting our /|^(4m^) and /|^(4m^) parameters to 

■^We use Eq. (^2|) to express I2''' in terms of I^^^ ^-nd Iq^ '' ■ 

''We remind to the reader that we have assumed a 5% error in the data of these two references. If 
we had assumed a 1(10)% errors, the statistical fluctuations quoted in Eq. ( |l3| ) would have decreased 
(increased) roughly by a factor of 5 (2). 
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Data Refs. 
J = / = 


i5|, ig 
J = 1 = 1 


Data Refs. 
J = / = 


0, |18| 
J = 1 = 1 




-0.0323 ±0.0005 - 


-0.1157± 0.0017 


-0.0289 ±0.0005 - 


-0.0980 ±0.0017 






0.011 ±0.006 




0.072 ±0.005 


/f (4m2) 


-0.031 ±0.004 




-0.003 ±0.003 






0.233 ±0.011 


0.0292 ±0.0005 


0.204 ±0.008 


0.0356 ± 0.0006 




0.237 ±0.007 


0.0050 ±0.0002 


0.267 ±0.006 


0.0062 ± 0.0002 


xVdof 
/2^''(4m2) 


1.6 

0.0082 ± 0.0008 


3.9 

-0.0016 ±0.0007 



Table 1: Fitted {I^{Am'^)) parameters to the experimental data of Refs. ||T|] ( J = / = 0), and Q 
[J = I = 1). J„(4m^) are given in units of (2m)" and /^'''(4m^) is calculated from Eq. (|l2). Errors in 
the fitted parameters are statistical and have been obtained by increasing the value of by one unit. We 
also give the threshold parameters a/j and hij obtained from an expansion of the scattering amplitude, 
ReT/,/ = — 167rm(s/4 — m?Y\aij + bij{s/A — m?) + • • •] close to threshold. In the the s— wave channel 
and due to the lack of error estimates in Ref . , we have assumed a rule of thumb error of 5% in the 
data and carry out the fit up to 900 MeV. We have chosen this set of data because the data of Refs. 
and seem to be inconsistent between each other at low energies and cover near threshold a narrower 
region than that covered by the analysis of Ref. For the p— wave, the data of Refs. jlSj (here again, 
we assume an error of 5% in the data, because the lack of error estimates of the original analysis ) and 
Ip^ disagree again, specially close to threshold. 



zero, and taking / and m as the physical parameters. Direct inspection of our expressions 
shows that one can not simultaneously absorb the a— and p— channel divergences by 
a redefinition of the m and / parameters. Since the model of Refs. [^-[|10] has been 



constrained not to have these new parameters (/|^4(4m^)) generated by the off-shellness it 
is impossible for them to simultaneously reproduce the a— and p— channels; a polynomial, 
with some new parameters, has to be added to the on shell-constrained solution of the 
BSE to properly describe the resonance region in the J = I = 1 channel. The inclusion 
of this polynomial with free parameters is justified @]-[|l3l within the lAM approach [^. 
We would like to stress here that considering explicitly the off-shellness automatically 
embodies the lAM, as it is invoked in Refs. and [|lOl, and thus incorporates such a 
polynomial. 



4 Conclusions and Outlook. 



We have solved the BSE for tttt scattering in the ladder approximation. This is to say 
using the potential and pion propagator at lowest order in the chiral expansion. This 
calculation produces unitary amplitudes in the s— channel and describes satisfactorily s— 
and p— wave phase-shifts from threshold up to the region of low-lying resonances. The 
present approach can be extended in principle to higher orders in the chiral expansion, i.e., 
including in the potential and in the pion propagator terms of order 0{p'^) and higher. 
The new divergences will become more severe and more subtraction constants will be 
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needed. That will be translated in an increasing number of free parameters, as it is the 
case also of standard ChPT. To be more specific, for instance, the 0{p'^) potential contains 
unitarity corrections in the t— and u— channels, (the corresponding s— channel correction 
to the amplitude is not two particle irreducible). Hence, the solution of the BSE with this 
potential will automatically implement exact unitarity in the s— channel and perturbative 
unitarity in the t— and m— channels. 

The inclusion of unitarity corrections in the t— and u— channels in the potential makes 
the practical solution of the BSE a cumbersome task and it is currently underway 



However, if one neglects these corrections in the t— and u— channels the BSE can be 
solved using a much simpler algebraic procedure as shown here for the lowest order case. 
In that case only unitarity in the s— channel will be restored and crossing symmetry will 
be violated. Thus, one should expect a general solution, below the four pion production 
threshold, of the form 

where P^"^(x) and Ql^{x) are polynomials^ of order n (corresponding to the order 0{p^^) 
in the proposed expansion) with real coefficients. Most of these coefficients have to be 
fitted to the data to accomplish with the renormalization program. The resemblance to 
a sort of pade approximant for the inverse amplitude, although not exactly in the form 
proposed in Ref. is striking. On the other hand, Eq. (|T^ provides a model of the type 
N/ D for the amplitude, where D has a right cut and the function N is approximated 



by a polynomial with no cuts. These ideas have been recently examined in Ref. P4 | 



Finally, we should mention that to describe the s— wave at higher energies than 0.8 
GeV, a coupled channel formalism needs to be used. Such an improvement of our model 
can be easily implemented, at least at the lowest order 0{p'^). 
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Appendix 

The 0{p^) + 0{p^) on-shell S'f/(2)— ChPT isoscalar s— and isovector p— wave tttt am- 
plitudes, expressed in terms of the renormalization invariant parameters /i,2,3,4; are given 

^Note that the first non-vanishing coefficient in Qn corresponds to the power (s — Am?)'-' . 
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by i 



TiAs) 



VjAs) 



Vl'{s)/f+{vl\s) + Us)x[Vl\s)] 

m?-2s 



2\ / /4 



(15) 



2 



I = J = 



I = J = 1 



1 1 , 



1 rf(cos^) 
a 2 



{fu{t)Ut) + /,j(n)/o(M)) Pj(cose) 



5m4 I imm'^{s-im?) . I9lm? {s-im?f . T — T — n 

1^ 96^ ^ 288;;^ ^ i - J - U 

89m^(s-4m^) _ 37m^(s-4m^)^ 



where t 

fiAx) 

9ij 



+ ■■• / = J = l 

-2(s/4 - m2)(l - cos 6*), m = -2(s/4 - m2)(l + cos 6*) and 



2887r2 28807r2 
.2\/i /3^ „, _ O/0//I ^2^ 



lOx^ + x{2s - 32m2) + 37m^ - Ssm^ I = J = 
2x'^ + x(s + 2m^) — m'^ — 4sm^ / = J = 1 

' (4011 + 8O/2 - ISlg + 8414 + 125) + m2(s - 4m2) x 
(3211 + 48r2 + 2Ah +'-f) + is- Am')' {fh + 1^2 + '-f) 



s—im 



Am' (-2/1 + 2/2 + 3/4 + 1) + (s - Am') (-2/i + 2/2 + 1 



(16) 

J = J = 
I = J=1 



On the other hand, expanding up to order (9(p^) the amphtudes of Eqs. (^) and (|T^) 
we find, 



Tij{s) = V,''{s)/f'+{W{s)+Io{s)x[V,''{s)] 
W'^(s) = 



2\ / f4 



(17) 



4^'^ + 7m'I^'' + Us- Am')I^'' + I^'' x [ V^oo(s) ]2 J = J = 



I = J=1 



Taylor expanding around threshold the function hjj and identifying the functions ¥2'^ 
of Eq. (|l^) and W/"' of Eq. (l7) one obtains the constraint given in Eq. (|1^) and the 
following relations: 



- 107 2/ll^rRp 2^^rRs 

^ 750 V 25 ° 25 ° 



(1^ 



1997 
3000 



7r 



2 f^^jR,p I 1152 ^^; g 

UOO ° 100 ° 
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- _ 2761 2 /1472 o,, 1776 o 224 o 64 

~ 3750 +'' U75^° 125^° + 75m2^2 + Sm^^^ j 



173 2 

120+^ VT^'' ' 3m2 
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